A matlab 1 package for exponential integrators
نویسندگان
چکیده
Recently, a great deal of attention has been focused on the construction of exponential integrators for semi-linear problems. In this paper we describe a matlab package which aims to facilitate the quick deployment and testing of exponential integrators, of Runge–Kutta, multistep and general linear type. A large number of integrators are included in this package along with several well known examples. The so-called φ-functions and their evaluation is crucial for stability and speed of exponential integrators, and the approach taken here is through a modification of the scaling and squaring technique; the most common approach used for computing the matrix exponential.
منابع مشابه
Norges Teknisk-naturvitenskapelige Universitet Expint — a Matlab 1 Package for Exponential Integrators
Recently, a great deal of attention has been focused on the construction of exponential integrators for semi-linear problems. In this paper we describe a matlab package which aims to facilitate the quick deployment and testing of exponential integrators, of Runge–Kutta, multistep and general linear type. A large number of integrators are included in this package along with several well known ex...
متن کاملExponential Rosenbrock integrators for pricing American options
In this talk, exponential Rosenbrock methods are investigated for applications in finance. The methods are shown to be an alternative to other existing procedures for the numerical valuation of American type contracts. Exponential Rosenbrock methods have good stability properties, they are fully explicit and do not require the numerical solution of linear systems as in contrast to standard inte...
متن کاملImplementation of exponential Rosenbrock-type integrators
In this paper, we present a variable step size implementation of exponential Rosenbrock-type methods of orders 2, 3 and 4. These integrators require the evaluation of exponential and related functions of the Jacobian matrix. To this aim, the Real Leja Points Method is used. It is shown that the properties of this method combine well with the particular requirements of Rosenbrock-type integrator...
متن کاملTime-stepping methods for the simulation of the self-assembly of nano-crystals in Matlab on a GPU
Partial differential equations describing the patterning of thin crystalline films are typically of fourth or sixth order, they are quasior semilinear and they are mostly defined on simple geometries such as rectangular domains. For the numerical simulation of these kind of problems spectral methods are an efficient approach. We apply several implicit-explicit schemes to one recently derived PD...
متن کاملA Krylov Subspace Algorithm for Evaluating the Φ-functions Appearing in Exponential Integrators
We develop an algorithm for computing the solution of a large system of linear ordinary differential equations (ODEs) with polynomial inhomogeneity. This is equivalent to computing the action of a certain matrix function on the vector representing the initial condition. The matrix function is a linear combination of the matrix exponential and other functions related to the exponential (the so-c...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2005